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ABSTRACT 

This  publication  describes  and  evaluates  a  computer 
program  useful  for  the  statistical  processing  of  bivar¬ 
iate  data  fields  expressed  in  polar  coordinates. 


FOREWORD 


An  IBM  7070  computer  program  it  described.  This  program  ie  used 
by  the  U.  S.  Naval  Oceanographic  Office  for  the  statistical  processing 
of  bivariate  data  fields  expressed  in  rose  form.  The  program  logic  is 
explained  and  the  quality  of  the  program  is  evaluated. 

Scientific  institutions  desiring  more  specific  information  concerning 
this  program  are  invited  to  write  to  the  U.  S.  Naval  Oceanographic 
Office. 


U.  S.  Naval  Oceanographic  Office 
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A  STATISTICAL  ROSE  PROGRAM 


I*  INTRODUCTION 

An  IBM  7070  computer  program  in  available  which  will  calculate 
five  parameters  of  frequency  distribution  over  a  two  dimensional  field 
(the  components  of  the  vector  mean,  the  components  of  standard 
deviation  along  the  principle  axes  of  the  distribution,  and  the  inclination 
of  the  major  axis),  given  either  raw  data,  or  vector  sums  and  vector 
sums  of  squares  of  data,  or  processed  roses.  Having  calculated  the  five 
parameters  of  distribution  density,  the  program  then  integrates  fre¬ 
quencies  in  the  polar  coordinates  of  the  bivariate,  tabulating  frequencies 
for  as  many  as  six  interval  groups  in  eight  cardinal  directions, 
assuming  normal  distribution. 


The  program  has  several  uses: 

« 

1.  It  may  be  used  for  the  statistical  processing  of  any  kind  of  rose, 
including  wind  roses,  sea  and  swell  roses,  current  roses,  cyclone 
roses,  etc.  Program  application  is  not  limited  to  velocity  fields.  The 
program  may  be  applied  to  any  bivariate  field  for  which  sampling 
distribution  may  be  normal  and  elliptical. 

2.  The  program  may  be  used  for  the  statistical  processing  of 
roses  for  which  data  are  insufficient  to  conclude  that  a  rose  based  on 
raw  data  is  representative.  The  program  calculates  the  parameters 
of  distribution  density  associated  with  the  data  rose  and  produces  a 
statistical  rose  with  the  same  parameters.  Assuming  that  the  data 
sampling  is  random  and  sufficient  for  calculating  representative 
parameters  of  distribution  and  that  an  infinite  random  data  sampling 
would  yield  a  normal  elliptical  frequency  distribution  field,  the  statisti¬ 
cal  rose  may  be  regarded  as  representative  for  the  rose  station.  In 
established  study  routines,  statistical  roses  have  been  generated 
assuming  that  samples  of  twenty  or  more  observations  per  station 
are  adequate  to  yield  representative  parameter  values  but  that  samples 
of  less  than  a  hundred  observations  would  not  be  adequate  to  yield  a 
representative  rose  using  standard  methods. 

3.  Under  certain  circumstances  the  program  may  actually  be  used 
to  calculate  statistical  roses  for  places  where  no  data  are  available. 
Since  the  program  is  designed  to  calculate  a  statistical  rose  from 
five  parameters  of  distribution  density,  these  five  parameters  may 
be  fed  to  the  program  as  data  to  yield  the  corresponding  statistical 
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roae.  For  an  area  with  scattered  stations  one  may  interpolate  by 
analysis  procedures  the  five  parameters  to  be  used  as  data  for  those 
places  with  no  data  sampling. 

4.  The  program  may  be  adjusted  to  yield  arbitrary  speed-interval 
classes  for  the  processed  rose.  The  data  rose  may  thus  be  processed 
for  speed-interval  classes  which  differ  from  those  given  with  the 
rose.  If,  for  example,  the  speed-interval  breakdown  of  Weather  Bureau 
data  differs  from  a  standard  speed-interval  format,  such  data  may  be 
processed  to  represent  standard  speed  intervals  by  employing  this 
program. 

5.  The  program  may  be  used  to  test  the  assumption  graphically 
that  large  bivariate  data  samplings  tend  toward  normalcy  of  distribution. 

II.  PROGRAM  SPECIFICATIONS 

The  program  was  written  in  the  IBM  7070  4-tape  autocoder  language 
so  that  maximum  use  could  be  made  of  efficient  logic  and  fixed  point 
arithmetic  operations.  The  program  occupies  3572  memory  storage 
addresses  at  the  time  of  this  publication,  including  the  IBM  input- 
output  subroutine  package,  and  data  fields. 

The  program  generally  requires  twelve  to  fifteen  seconds  to  process 
a  rose.  The  frequency  integration  routine  accounts  for  most  of  the 
time  required  for  the  production  of  a  rose.  Only  a  few  milliseconds 
are  required,  in  general,  for  calculating  the  distribution  density  para¬ 
meters  of  a  rose,  assuming  that  the  data  rose  input  format  is  used. 
Since  the  machine  rents  for  $2.50  per  minute,  machine  cost  for  a 
single  rose  ranges  between  fifty  and  seventy  cents. 

Special  routines  written  for  the  IBM  1401  are  used  for  reading  data 
cards  and  listing  the  completed  rose. 

in.  PROGRAM  LOGIC 

A.  Calculating  the  vector  mean,  vector  variance  and  covariance: 

1.  “Raw"  data:  Given  raw  data  consisting  of  two  vector  compo¬ 
nents  per  sample  (cartesian  or  polar  coordinates,)  a  subroutine  may 
be  employed  which  calculates  the  sum  of  components,  the  sum  of  the 
squares  of  components  and  the  sum  of  component  cross-products  in 
the  cartesian  coordinate  system. 
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Suppose  that  a  raw  data  format  is  available  which  expresses  wind 
data  in  degrees  from  north  and  knots.  The  subroutine  sums  and  counts 
successive  data  samples  and  calculates  the  components  of  the  vector 
mean,  vector  variance,  and  covariance  according  to  the  formulae: 


v  =  v.  sin  9 

xj  j  j 

v  s  v,  cob  9, 

yj  j  j 


(1) 

(2) 


7J  p 


xj 


VY  *  N  IVyj 


S  2  a  £v  “  -  V 


2  .  y  2 

X  N  L  x j  X 


e  2  1  2  -  2 

*  —  tv  -  v 
y  N  2.  yj  y 


cov  s  i  rv  v  -  vv 


N  IVxj 


yj  x  y 


(3) 

(4) 

(5) 

(6) 

(7) 


where 


N:  Total  frequency  of  observations 
j  :  The  jth  data  sample 
9^:  Degrees  from  north  of  the  jth  data  sample 

VjS  speed  (knots)  of  the  jth  data  sample 


x:  The  easUwest  coordinate  suds 
y:  The  north-south  coordinate  axis 


vx:  The  X- component  of  the  vector  mean 
v  :  The  Y- component  of  the  vector  mean 

y 
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Sx  :  The  X-  component  of  the  vector  variance 
2 

Sy  :  The  Y- component  of  the  vector  variance 
coy :  The  covariance  of  the  X-Y  components 


2.  Data  rose:  Given  a  data  rose,  wherein  raw  data  has  previously 
been  tabulated  in  specified  speecUinterval  classes  for  the  eight  cardinal 
directions:  the  vector  mean,  vector  variance  and  covariance  may 
be  calculated  using  the  formulae: 


r*  -  x\  Fw 

IVd2  *  Iai2Fid  (d  «  N,  NE,  E.....NW) 

;,4{|rE-iv,0-10,u 

['IVNE  -  Xvsw,+  (  IVSE  -  rVNw']} 

;y'Tf{'ZvN-Ivs>t  0-707U 

[<XVNE  -  XYSW»  -  (XVSE  -  XVNw']} 

VV  ■fjpX’E2  *  Xyw2>  +  °-5  [ 

5t2ne  *  Xy2sw» +  <Xy2sE  +  Ey2W]} 

Vy  '  M'Xvn2  ♦  Ivs2'  +  °-5  [ 

«Iy2NE  *  Xy2SW»  +  <Xy2SE  *  Xy2NW»]  } 

Vy  •  TT  {  0,5  [(Xy2NE  +  Iv2sw>  -  (Iv2se  +  Iv2nw>]}<9» 


then 


.  a*  .  6 

S  ■  v  v  -  V 

x  xx  x 

_  2  -  2 

S  *  v  v  -  v 

y  y  y  y 


COV  ■  V  V  -  V  V 

X  y  X  y 


(10) 
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where  i :  the  ith  speed  group  in 
d:  the  dth  direction  class 

a,  :  the  midpoint  speed  in  the  ith  speed  group 

F.^ :  the  frequency  of  observations  in  the  (i-d)th 
speed-direction  group 


Because  the  centroid  of  distribution  in  a  speed-direction  group 
is  not  generally  coincident  with  the  point  of  intersection  of  the  group 
midpoint  speed  and  the  ray  bisecting  the  octal  direction  class,  this 
method  of  calculating  distribution  parameters  may  not  yield  wholly 
accurate  results.  This  method  characteristically  exaggerates  the  distri¬ 
bution  spread,  for  instance,  producing  an  overestimate  for  the  scalar 
mean  velocity  and  the  vector  variance.  Functions  may  be  applied  which 
will  approximately  correct  for  the  interval  errors  intrinsic  in  this 
method;  however,  the  program  does  not  presently  incorporate  these 
correction  function*. 


B.  Calculating  the  inclination  of  the  principle  axes  of  the  distribution 
and  the  components  of  standard  deviation  along  the  principle  axes 
of  the  distribution:  The  equations  used  for  calculating  the  parameters 
of  the  distribution  were  derived  from  equations  provided  in  “An 
Introduction  to  the  Theory  of  Statistics."* 


Let 

_  2  _  2 
=  S  +  s 
x  y 

(ID 

K2 

e  2  _  2 
=  S  -  S 
x  y 

(12) 

L 

2  2  "i- 

a  (K2  +  4  COV  )  L 

(13) 

then  E^ 

«  [o.5(K1  +  L)]“^ 

(14) 

E2  *  (K1  *  E1  '  <15) 


*  G.  U.  Yule  and  M.  G.  Kendall,  “An  Introduction  to  the  Theory  of 
Statistics,”  Edition  13;  Charles  Griffin  fc  Company,  Limited;  London, 
1949*  chapter  12,  pp.  227-232. 
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Where 


E  :  The  component  of  standard  deviation  along 
the  major  axis  of  the  distribution 

E^  :  The  component  of  standard  deviation  along 
the  minor  axis  of  the  distribution 

Let  ••  represent  the  angle  included  between  the  X-coordinate  and  one 
of  the  principal  axes  of  the  distribution.  Then 

tan  28*  a  2cov  / 

or  8*  a  0.5  arctan  (2cov/ K^)  (16). 


If  >0,  let  8  a  8*  ;  otherwise  let  8  a  8*  +  90*.  If  K  a  0,  then  if 
cov<  0,  8  a  .45*;  otherwise  8  a  +45*.  Then  8  identifies  tne  major  axis 
of  the  distribution. 


C.  Integrating  frequencies  for  the  normal  surface  within  the 
limits  of  specified  direction  and  speed  classes 


1.  General  description  of  the  method  of  integration:  The  program  . 
regards  the  normal  surface  as  a  grid  of  2,228  areas  of  known  frequency. 
The  grid  is  uniquely  defined  by  the  parameters  of  distribution  density. 
The  program  locates  the  centroid  of  each  such  area  in  the  polar 
coordinates  of  the  bivariate  and  accordingly  tabulates  the  known 
frequency  in  the  appropriate  speed-direction  group. 


2.  Assumptions  employed  in  deriving  the  grid  of  areas  of 
known  frequency:  The  normal  surface  may  be  visualized  as  a  bell 
shaped  mound  with  elliptical  horizontal  cross-sections,  its  axis  of 
symmetry  representing  the  ordinate  which  measures  density  of  distri¬ 
bution  frequency  and  perpendicular  to  the  plane  which  measures  vector 
variates.  The  equation  for  distribution  density  (the  normal  surface) 
is  given*  for  a  point  (x,y)  in  cartesian  coordinates  coincident  with  the 
principal  axes  of  the  distribution: 


+  (y/®2)2 

x.y  *  ,vr  F  - 


2TElE2 


(17.a) 


*  Yule  and  Kendall, op.  cit. 
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The  exponent  in  equation  (17*a)  is  elliptical  in  form.  Let  a  cross* 
section  of  the  distribution  be  identified  with  the  parameter,^,  such 
that 


(18.a) 


The  parameter  ^  identifies  an  ellipse  which  intersects  the  major  axis 
of  the  distribution  at  the  point  Ey  and  the  minor  axis  at  the  point 
E^.  X,  then,  is  a  factor  of  proportionality  relating  elliptical  cross* 
sections  of  the  distribution,  and  is  everywhere  constant  along  the 
perimeter  of  an  ellipse.  It  follows  from  equation  (17.a)  that  the  density 
of  distribution  is  likewise  constant  at  each  point  along  the  perimeter 
of  an  elliptical  cross. section,  and  that  it  may  be  given: 


.4* 

21TE  E 
12 


(17.b) 


Let  the  polar  coordinates  (R,0)  be  defined  such  that-£  =  cos0;X  * 
sin0.  Then  x  a  Rcos0;  y  ■  R  sin0,  and  the  ellipse,  equation  (lf.a) 
may  be  expressed  in  the  polar  coordinates  (R,0): 

2  2  2y2  ,  2  2 ,  _  2  2-  , 2  y 2  iiq 

R  a  Ej  E2  f  /  (E2  cos  0  +  Ej  sin  0)  a  IE  \  (18.b) 

and  x  >  IE  |  cos0,  y  a  IE  ^  sin0 
The  frequency  of  distribution  is  given  for  a  segment: 


It  is  convenient  to  integrate  equation  (19. a)  between  ellipses.  Then 
the  integral  (19.a)  may  be  transformed  to  the  elliptical  coordinates, 
(*.*>: 
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*\*m 

0  * 

d0 

(19.b) 

/x  y\ 

where  the  Jacobian,  J\>  0  )  =  <J  x  dy  - 

dx  i 

ly  =  IE2^ 

(20) 

if  ii 

00  "i 

•  0  that  Fj.  ,  »  f  E2  )[ 

d0 

(19.c) 

Since  ^  and  0  are  mutually  independent,  functions  involving  ^  only 
may  be  extracted  from  equation  (19.c)  as  constants  with  respect  to 
0,  and  conversely.  Then,  integrating  between  the  limits  K  and  \  , 

*,  “<■  *r 


Evaluating  equation  (22)  between  the  limits  =  0,  ^  31  \)  aad 
(0j  »  0,  02  =  2 ir)  provides  the  frequency  of  distribution  enclosed 

by  the  ellipses: 


F(*)  =  1  -  e  2 


(23) 


Conversely,  given  the  frequency  of  distribution,  F,  enclosed  by  an 
ellipse,  to  identify  thej^  parameter  of  the  ellipse: 

*  =  [-2  ln(l  -  F)]7  (24) 

It  may  be  noted  that  the  0  integral  in  equation  (21)  provides  the 
area  of  the  sector  of  the  unit  ellipse  which  is  subtended  by  the  angle 
(02  -  0}).  The  ^  integral  provides  the  frequency  of  distribution  between 
the  ellipses  ^ _  and  X  .  The  equation  (21)  therefore  shows  that  sectors 

A  mt 
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including  equal  areas  of  the  unit  ellipse  likewise  include  equal  fre¬ 
quencies  of  the  distribution  in  the  ring  between  adjacent  ellipses.  The 
normal  elliptical  distribution  may  therefore  be  subdivided  into  seg¬ 
ments  of  known  frequency,  employing  a  grid  of  sectors  subtending 
equal  areas  of  the  unit  ellipse  and  ellipses  enclosing  known  frequencies 
of  the  distribution. 

The  grid  of  sectors  of  equal  area  may  be  derived  by  employing  an 
arbitrary  parametric  angle,  ^  ,  such  that  points  on  the  unit  ellipse 
may  be  given  in  cartesian  coordinates  coincident  with  the  principal 
axes  of  the  distribution: 


accordingly, 


P  s  (Ej  costfr,  sin^f) 
E2 

tan$»=r 


(25) 


so  that  the  area  of  the  sector  defined  by  is  given: 

A  =  flf?  arctan  (  fl_  tan  t)  *  flfz  * 

2  '  E£  '  2 

and  integer  multiples  of  the  parametric  angle  ^  shall  define  integer 
multiples  of  the  unit  sector  area  defined  by  \ff  . 

Given  the  grid  of  segments  of  known  frequency,  it  is  required  to 
find  a  representative  point  in  each  segment  such  that  if  the  segment 
overlaps  adjacent  speed-direction  groups  the  point  shall  be  tallied  in 
the  group  containing  the  greatest  portion  of  the  segment  frequency. 
The  centroid  of  the  segment  may  be  regarded  as  the  desired  repre¬ 
sentative  point.  The  coordinates  of  the  centroid  of  a  segment  may  be 
determined,  in  cartesian  coordinates,  (x,y),  coincident  with  the  principal 
axes  of  the  distribution: 
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Since  x  s  coi  ,  y  =  E2  k  sin  ^ 

the  equations  (26)  may  be  transformed  to  the  elliptical  coordinates, 

(f.  tb 


T  */>  /E i  *cos*  r*  J(t  *) d*  d* 

d* 

Ei  */*/'*  **  c°8 *  d*  d* 

"  ■  l/lf-n'-r* 


•  E>  fib 


n  \jf  J cot  \jf  d  jr 


/d^f 


JL/ 

H(J)]  =  j/e  2  |2  dj,  F(I)J  =  j/< 

Likewise,  f  =  E  H(  i)] 


where 


-J-*2 

:2'  *d* 


*71)  J, 


(27) 


(28) 


F(^),  evaluated  between  ellipses  *“d  %  ,  provides  the  frequency  of 
distribution  between  the  two  ellipses.  The  function,  H( 4^),  may  be 
evaluated  between  the  same  limits: 


2n+3 


2n+3 
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) 


2n  (2n+3)  n J 


(29) 
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Let  (n  -  2^1  define  successive  rays  including  sectors  of  equal  areas 
in  the  unit  ellipse  for  n  =  0,  1,  2,  etc.  Then  integrating  the  equations 
(27)  and  (28)  between  the  limits  ({j.l)>  (^j)and  (n-£)[j/,  (n+^)\f/  provides 
the  components,  Pj,n,  of  the  centroid  of  the  corresponding  segment: 


sin(n-|)^  E2  [cos(n-j)^ 


*  (Ej  cos  n\ft,  E2 

where  7j  *  £|||J  ^  G(^).  and  G(^)  =  2  sinj^ 


(30) 


*j-l 


That  is,  the  ray  defined  by  n^r  bisects  the  area  of  the  sector  defined 
between  the  rays  (n-4)  ^  and  (n+£-)  \fk  and  so  bisects  the  frequency  of 
distribution  included  in  the  sector.  It  may  be  noted  that  when  n  *  0, 
P  Q  i_a  proportional  to  the  module,  Ej,  of  the  major  axis.  When  n(|rs 
90»,  P.  is  proportional  to  the  module,  Ei  of  the  minor  axis.  The 
locus  o/P.  for  j,^ constants,  is  an  ellipse. 

The  function,  G(^)  is  the  factor  providing  correction  for  sector 
parallax.  It  may  be  noted  that  as  the  parametric  angle,  approaches 
the  limit,  0,  G(^)  approaches  the  limit,  1,  so  that  the  sector  parallax 
error  vanishes  as  the  sector  angle  vanishes. 

The  parameters  ^  j  and  n ^  thus  define  the  grid  of  segment  centroids, 
and  may  be  predetermined.  The  equation  (30)  relates  the  grid  to  the 
individual  rose  operated  upon  by  the  program.  Grid  choice  is  arbitrary. 
Grid  segments  were  defined  for  this  program  so  as  to  enclose  discrete 
values  of  distribution  frequency.  The  corresponding  {  j  values  were 
then  calculated  employing  equations  (25),  for  given  n- intervals,  and 
(24).  The  grid  was  chosen  for  segment  definition  such  that  segment 
areas  would  be  of  comparable  magnitude  and  of  similar  shape,  and 
such  that  maximum  programming  advantage  could  be  taken  of  grid 
symmetry.  The  grid  specifications  are  presented  in  table  I. 
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3.  Locating  the  segment  of  known  frequency  in  the  grid  of  direction- 
speed  interval  classes:  The  equation  (30)  expresses  the  point,  Rj  n, 
in  cartesian  coordinates  coincident  with  the  principal  axes  of  the 
distribution,  and  is  transformed  to  components  in  the  coordinate 
system  of  the  bivariate  with  the  equations: 


Xj,n  =  (E^osO)  ( cos  nity  -  (E2  sinO)  (7.  »in  n^r)  +  vx 

Yj  n  *  (E^cosd)  ( ^  sin  nf)  +  (E^inO)  (T.cos  nty  +  vy  (31) 


The  point  is  located  within  the  limits,  Vj  and  V2,of  a  designated 
speed  interval  class  if 


<  (X 


j.n 


+  *> 
J»n 


<  V, 


The  point  is  located  within  the  appropriate  octal  direction  class 
according  to  the  following  logic: 


located  octant 


(1)  If 

IVl*“22i‘  >b.» 

,  then: 

(a)  if  X  >0. 
J»n 

EAST 

(b)  otherwise. 

WEST 

(2)  M 

>|XJ.n 

|  ,  then: 

(a)  if  Y.  >0, 
J.n 

NORTH 

(b)  otherwise. 

SOUTH 

(3)  Otherwise,  if  X.  >0, 

J*n 

then: 

(a)  if  Y  >0. 

NORTHEAST 

(b)  otherwise. 

SOUTHEAST 

(4)  Otherwise,  X.  <  0, 

J*n 

then: 

(a)  if  Y,  >0, 

NORTHWEST 

(b)  otherwise. 

SOUTHWEST 
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IV.  ROSE  QUALITY  EVALUATION 


It  was  not  practical  to  hand-calculate  a  rose  for  control  results 
after  the  manner  employed  by  the  program.  Hand-calculation  of  a 
single  rose  would  have  required  several  months  to  perform,  and  would 
at  best  have  provided  a.  check  against  the  arithmetic  of  the  program 
rather  than  a  check  against  fundamental  program  logic.  Rather,  test 
data  known  or  assumed  to  be  normally  distributed  were  processed  by 
the  program  and  the  results  were  examined  for  correspondence  with 
the  test  data. 

A.  Brooks'  circular  normal  rose  data 

Brooks'*  roses  supplied  from  tables  were  used  as  one  source 
of  test  data.  These  roses  are  theoretically  characterized  by  normal- 
circular  distribution  and  are  useful  for  testing  the  program  logic, 
equation  (30),  for  the  special  case,  =  E £.  The  table  roses,  however, 
are  derived  by  graphical  interpolation  methods,  and  may  not  be  exactly 
normal  nor  circular  in  fact.  A  representative  test  sample  is  provided 
in  appendix  I.  The  Brooks'  data  rose  is  listed  along  with  the  components 
of  tiie  vector  and  scalar  means  which  theoretically  identify  it.  The 
program  image  of  the  data  rose  is  listed  for  easy  comparison. 

B.  Program  roses  as  test  data 

Theoretically,  the  program  supplies  roses  which  are  elliptically 
and  normally  distributed.  Consequently  the  program  rose  supplies  data 
ideal  for  a  test  against  itself.  A  representative  sample  test  is  provided 
in  appendix  n. 

C.  Analysis  of  test  results:  This  test  provides  an  index  of  rose 
degeneration  due  both  to  the  method  of  rose  frequency  integration  and 
to  the  method  (Section  IU.A.2)  for  calculating  the  parameters  of  the 
data  rose.  The  method  of  rose  frequency  integration  may  be  expected 
to  produce  distortions  in  the  image  rose  because  of  segment  area 
overlap  between  adjacent  speed-direction  groups.  However,  the  fre¬ 
quency  tabulation  error  due  to  the  overlapping  of  a  segment  between 
adjacent  speed-direction  groups  may  be  either  positive  or  negative 
for  a  group,  and  the  collection  of  such  errors  tend  in  general  to  cancel 


*  Brooks  and  Carruthers,  "Handbook  of  Statistical  Methods  in  Meteorol¬ 
ogy,"  M.  O.  538.  London,  Her  Majesty's  Stationery  Office,  1953; 
Chapter  11. 
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as  the  frequency  population  of  the  group  is  tallied.  The  method  of 
frequency  integration  therefore  probably  doee  not  account  for  conaiatent 
divergencea  between  the  data  roae  and  its  program  image.  Diatortiona 
due  to  thie  method  may  be  likened  to  a  “white  noise"  which  ia  dimini ahed 
in  intenaity  aa  the  frequency  grid  ia  divided  into  emaller  and  amaller 
aegmenta.  The  method  (Section  UI.A.2),  on  the  other  hand,  may  be 
expected  to  produce  image  roae  valuea  for  the  acalar  mean  and  vector 
deviations,  and  £2*  which  conaiatently  exceed  the  correaponding 
valuea  for  the  data  roae. 

In  appendix  I  it  may  be  noted  that  £}  and  £2  are  only  approximately 
equivalent.  Thia  may  be  taken  to  indicate  either  that  the  object  roae 
does  not  represent  an  exactly  circular  distribution  or  that  the  method 
for  calculating  the  parameters  of  the  data  roae  may  thus  distort  the 
relationship  between  theae  parameters.  The  more  probable  alternative 
is  that  the  object  roae  does  not  represent  an  exactly  circular  distribu¬ 
tion;  in  appendix  11,  wherein  a  program  rose  ia  compared  with  its  own 
image,  it  may  be  noted  that  the  inclination  of  the  major  axis  and  the 
ratio  of  deviations  between  the  principal  axes  closely  coincide. 

The  distortion  due  to  the  method  of  calculating  parameters  of  the 
data  rose  is,  however,  apparent  in  both  the  Brooks'  and  program  rose 
samples;  that  is,  the  program  image  values  for  the  scalar  mean  velocity 
and  the  vector  deviations  exceed  the  correaponding  parameters  for  the 
test  roae.  Nevertheless,  corresponding  values  between  the  test  rosea 
and  their  program  images  are  clearly  similar,  particularly  for  the 
program  roae.  Such  similarity  suggests  that  the  program  logic  is 
successful,  and  that  rose  distortion  due  to  segment  overlap  error  is 
small.  Unfortunately,  the  precise  extent  of  segment  overlap  error  is 
masked  by  the  distortion  which  may  be  attributed  to  the  method  for 
calculating  data  rose  parameters.  The  error  introduced  by  the  program 
method  for  frequency  integration  cannot  be  estimated  with  any  precision 
from  these  test  rssults. 

D.  Conclusion:  The  program  can  probably  be  improved  by  applying 
an  interval  correction  function  to  the  method  (Section  in.A.2)  for 
calculating  the  parameters  of  a  data  rose.  Doubtless,  too,  there  is  a 
more  efficient  method  for  integrating  frequencies  of  the  distribution  in 
the  speed-direction  interval  classes.  However,  test  samples  suggest 
that  the  program  rose  is  of  substantial  quality,  and  that  the  program 
may  be  employed  confidently  whenever  program  application  is  war¬ 
ranted. 
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Appendix  Is  BROdCS8  TSPT  DATA  ROSE 
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Appendix  II*  PR0Q8AM  TEST  DATA  BOSS 
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